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ABSTRACT 


The single terminal and two terminal diversity precipitation 
attenuation distributions are calculated for a millimeter wavelength 
earth -satellite propagation path. The calculated probability 
distributions are compared to attenuation probability distributions 
measured at two locations with the ATS-5 Millimeter Wavelength 
Experiment 15.63 GHz downlink. Using an empirically derived three 
mode exponential cell diameter - rain rate function, the calculated 
distributions agreed within -2 to +1 dB of the measured distributions 
for both terminal locations. The sensitivity of the calculated distri- 
bution to variations of the cell model parameters was also investigated. 
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CHAPTER I 
INTRODUCTION 


The need for more satellite communications channels has 
recently stimulated interest in the millimeter wavelength portion 
of the spectrum. The approval of new channels above 10 GHz by the 
World Administrative Radio Conference[l] and the proposal of several 
domestic satellite communications systems[2] in 1971 indicate the 
extent of this interest. However, before such frequencies can be 
used efficiently, more information is needed about the propagation 
phenomena which occur at these wavelengths. It has been known for 
some time that atmospheric components can attenuate millimeter 
wavelength signals significantly, as shown in Fig. 1. The greatest 
attenuation is caused by high liquid precipitation rates occurring 
along the propagation path. However, there is little data available 
about the effects of precipitation which are useful to an engineer 
designing a millimeter wavelength earth-satellite communications 
link. To study these precipitation effects and develop propagation 
prediction techniques for system design, a millimeter wavelength 
experiment which included path attenuation measurements at 15.3 GHz 
and 31 .65 GHz was included on-board the Applications Technology 
Satellite, ATS-5[3] . 
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Fig. 1 .--Attenuation due to gas absorption, fog and mist, and 
rain at millimeter wave frequencies (from Mondre[26]). 
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The attenuation probability distribution function is frequently 
used to predict long term effects of precipitation. This function 
gives the percentage of time that a particular level of attenuation 
is equalled or exceeded or, conversely, the system margin required 
to insure that the system is operational for a given percentage of 
time. Thus, it is useful in determining the economic feasibility 
and link reliability in system applications. This distribution 
function is strongly dependent upon the climatology of the region in 
which the earth terminal is located. Thus, the system designer needs 
a technique for predicting the precipitation attenuation probability 
distribution function for various earth terminal locations. 

A simple, mathematically tractable storm cell model is used 
to make precipitation attenuation distribution predictions for a 
single earth terminal and a two-terminal diversity configuration. 

The calculated distributions are integrals of a rain rate probability 
density and an intersected path length probability distribution. The 
rain rate probability density is determined from tipping bucket rain 
gauge measurements. Thus, the statistical characteristics of pre- 
cipitation events in different locations are included in the calculated 
attenuation distribution. A technique for using National Weather 
Service rain measurements is described so that easily available data 
can be used. The intersected path length probability distribution 
is calculated using a right circular cylinder with homogeneous rain 
rate, variable diameter and fixed height as a storm cell model. In 
this way, as a storm cell model, the physical characteristics of the 
precipitation events are represented. 
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Previous models used to estimate attenuation caused by rain 
have been based on postulating or measuring a precipitation distri- 
bution in a vertical slice along the propagation path. Dutton[4] 
assumed a spherical storm cell with a Gaussian variation of precipi- 
tation with altitude. The resulting attenuation calculated for a 
single fixed cell was an integral over the propagation path and a 
function of cell center location. This model is of limited use 
because it is difficult to incorporate the features which are needed 
to calculate the attenuation distribution. First, the statistical 
characteristics of the rain rate which must be added are hard to 
include. Second, the ways in which the cell might intersect the 
path to produce a particular attenuation are difficult to determine 
because of the vertical inhomogenity of the rain rate. Thus, a 
simpler model is required. A model proposed by Bradley[5] postulated 
a precipitation rate versus time model for an average storm in the 
geographic area of his interest. This model was then used to esti- 
mate the mean time between peak attenuation occurrences for a 
horizontal path. Although this approach uses the precipitation 
characteristics of a particular location, application is limited 
to horizontal paths. 

The approach which is used here is to lump the effects of all 
liquid precipitation events together and model them with a single 
cylindrical storm cell. This model is described and the attenuation 
probability distribution function for a single terminal is calculated 
in Chapter II. In Chapter III, the model is extended to a two 
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terminal spatial diversity configuration and the joint attenuation 
probability distribution is calculated. The parameters which are 
used to relate the cell model to the physical world are discussed 
in Chapter IV. In Chapter V, calculated results are compared to 
propagation data measured by two terminals using the ATS-5 15.63 GHz 
downlink. The sensitivity of the model to parameter variations is 
also discussed and a diameter - rain rate function is proposed which 
improves the agreement of the calculated results. 
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CHAPTER II 

CALCULATION OF THE SINGLE TERMINAL ATTENUATION 
PROBABILITY DISTRIBUTION 


As described earlier, previous methods for predicting attenu- 
ation due to rain have all used some type of rainstorm model. But 
the models proposed could not be readily applied to solving the 
communications design problem because they were too specific and 
could not be used in other geographic regions. Another practical 
constraint which is imposed by the design problem is that the model 
should be mathematically tractable. 

A reasonable approximation for the shape of a rainstorm is a 
cylinder. Thus, a right circular cylinder is used to model the 
precipitation phenomenon. To maintain the simplicity of the model, 
the precipitation in the cell is assumed to be in the form of rain 
only and the precipitation rate, r, is assumed to be non-zero and 
uniform throughout the cell. The cell diameter, denoted by d(r), 
is a function of the rain rate and represents an effective cell 
diameter. Consequently, the storm cell diameter may not compare 
well with observed rainstorm cell dimensions. 

The origin of the coordinate system will coincide with the 
earth terminal location on a flat earth such that the terminal- 
satellite propagation path lies in the x-z plane. The elevation angle. 


6 



will be defined from the x-axis toward the z-axis. The proposed 
rainstorm model with the coordinate axes is shown in Fig. 2. The 



SATELLITE 



Fig. 2. --Earth terminal and storm cell model. 


center of a storm cell will be defined as the point (x 0 >Y 0 ) where 
the vertical axis of the cell intersects the x-y plane. The 
occurrence of a storm cell centered at a point within the universe 
of storm cell centers is assumed to space-time ergodic. Thus, the 
time statistics of cell center occurrences observed at a point are 
equivalent to the spatial statistics of cell center occurrences at 
a particular point observed over the universe of cell center locations. 
The probability of a model cell centered at a point is assumed to 
be uniform over the universe of cell centers. 

Before calculating the single terminal attenuation probability 
distribution function, the physical processes involved will be 
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considered to determine the variables whose probability functions are 
needed. For non-zero attenuation to occur along the propagation path 
a section of the propagation path must be intersected by the storm 
cell. The rain rate along this intersected section of the propagation 
path and the length of the intersected path, i, determine the magni- 
tude of the attenuation according to the Gunn-East equation[7]. The 
intersected path length can be calculated directly from the inter- 
section geometry. Thus, to compute the attenuation probability 
distribution function, the storm cell rain rate probability density 
function and an intersected path length probability density function 
are needed. 

Therefore, for the random variables r and i, define 


( 1 ) 

and 


P r (r) = probability density function of the rain rate 
r in a model storm cell 


(2) p U) = probability density function of the length of 

the propagation path intersected by the storm 
cell. 

Since p r (r) and P £ U) are probability density functions then 


(3) 

and 


oo 

P r (r) dr = 1 
o 
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(4) 


£ max 


J 

o 


P z (i) d z = 1 , 


where z „ is the maximum intersected path length, 
max r 

The Gunn-East relation[7] mentioned previously gives a 
functional relation between the random variables, r and z, and the 
attenuation, a, in dB, 

(5) a = kr P £. 

The constants k and p are frequency and temperature dependent. 

Using a functional transformation theorem[8] for functions of two 
random variables it is found that for the two functions 

(6) a = kr P £ 

(7) r = r 
that 

(8) P a , r ( a ’ r) = P r#1 (r.O/|J(r.O| 

where p (a,r) is the joint probability density function of a rain 
ct,r 

rate r occurring on the intersected portion of the propagation path 
and the total attenuation occurring along the propagation path is a. 

The function p 0 (r,£) is the joint probability density function that a 
rain rate r occurs on the intersected path and the length of the 
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intersected propagation path on which the rain rate occurs is z. 
Finally, J(r,n) is the Jacobian of the functions a and r. The 
intersected path length, z, shown in Eq. (8) is the solution of 
Eqs. (6) and (7) in terms of a and r, 

(9) z = a/kr p . 

From an examination of the physical processes which are modelled by 
the storm cell it is clear that the occurrence of the rain rate, r, 
is independent of the occurrence of the intersected path length, z. 
Thus , 

00) P r>£ (r,£) = p r^ 

The Jacobian in Eq. (8) is given by 

01) J(r,£) = -kr p . 

Substituting Eqs. (9 )- (11 ) into Eq. (8) the result is 

02) P a ,r (a ’ r) = p r^ Pfc(aAr p )/kr p . 

If a storm cell having finite dimensions is considered, there 
is a maximum path intersection length, £ max ; and, consequently, for 
a given attenuation, a, there is a minimum rain rate, r min (a), which 
is 

(13) r m . n (a) ■ (»/ ktmax ,V p . 
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Thus, the attenuation probability density function, p a (a), can be 
found from p^ r (a,r) by integrating r over the range r^ n (a) to « 
which yields 


(14) 


CO 

p 0 (a) = [ p f (r) pJa/krPj/krP dr. 

* t \ 


r • _ (a) 
mi n 


The attenuation probability distribution function P a (<* 0 ), which is 
defined as the probability that the attenuation, a Q , is equalled or 
exceeded is given in Eq. (15). 


(15) 


W ■ 


P a (a) da. 


Substituting then p a (a 0 ) is 


CO oo 


(16) 


p (a.) 
a 0 


= 1 


a n r'Aa) 
o mi n 


P r ( r ) P a (a/kr p )/kr p dr da. 


To obtain an expression for p a (a 0 ) in terms of r and z only let 


(17a) 

(17b) 

(18a) 

(18b) 


x 

y 

u 

V 


r 

kr p ji 

r 

z 


and use Eq. (19) for transforming variables in a multiple integration 
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(19) 




f(x(u,v),y(u,v))|J(x,y)| dudv. 


where J(x,y) is the Jacobian, R is the region of the xy integration 
and R' is the corresponding region of the uv integration. Performing 
the indicated operations the result is 


max 


( 20 ) 


P o (o 0 ) = J j P r (r) P £ U) d£ dr 

r m -,„(0 *-i (<*„>!'') 
mi n o l o 


with the new function £-](a Q ,r) defined as 


(21) .£ (a ,r) = minimum intersected path length needed to 

produce an attenuation a Q given the rain 
rate, r. 


The function 2 ,^(a Q ,r) is a consequence of transforming the inte- 
gration from the a,r plane to the r,£ plane and is given as 

(22) MV r ) = a 0 /kr P 

and r miW is given by 

(23) r • (a ) = (a /k£ )^P. 

imn o o max' 

The £ integration defines the probability distribution function, 

P £ U), which is the probability that the intersected path length 

equals or exceeds i. Thus, the final form for P (a ) is 

a o 
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(24) 


P '(a) 
a 0 


= N 


r mir/ a o^ 


P r (r) P^(£ 1 (a o ,r)) dr. 


where N is a normalization constant which is included to insure 
that P a (o) = 1. Thus, it is given by 


(25) 




P r (r) p £ (*i (0.D) dr. 


In order to use Eq. (24) for numeric calculations, two 
probability functions, p r (r) and P^U), must be calculated. It is 
clear from the definition of P^U) that it is calculated from the 
storm intersection geometry. The rain rate probability function 
P r (r), will be determined from meteorological measurements made in 
the geographical region of interest. 

One way in which rain rate probability functions are measured 
is through the use of a tipping bucket rain gauge. Briefly, this 
gauge measures the time that it takes for a fixed amount of rain to 
fall. The number of tips occurring in a given time period then 
determine the average rain rate over that time period. The National 
Weather Service generally uses gauges which record a tip when 0.01 
inch of rain falls and the number of tips occurring in a sixty 
minute period beginning on-the-hour are recorded and publ ished[9]. 
Since these measurements are available for weather stations over 
the entire United States, the tipping bucket rain rate probability 
density function p B (r), which is defined as the probability density 
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of the rain rate r measured by a tipping bucket, will be used to 
calculate p r (r). In this way the geographic variation of precipi- 
tation phenomena is included in the attenuation probability 
calculation . 

Before this is done, it should be observed that the sixty 
minute averaging time used in the National Weather Service measure- 
ments is too coarse for accurate estimations of P a (a ). However, 
a method will be presented later which permits the estimation of a 
clock-minute rain rate probability density function from the clock- 
hourly rain rate probability distribution function. Although, it 
would be best to have the instantaneous rain rate probability density 
function, currently the clock-minute function is the finest scale 
that is available. 

If the measuring aperture of the tipping bucket rain gauge 
is small compared to the smallest spatial variation of the rain rate 
within a rainstorm, then, to a good approximation, the density 
function measured by the tipping bucket is a point function. In 
terms of the model storm cell it is seen that a non-zero rain rate 
is measured only when the center of a storm cell occurs within a 
circle with the same diameter as the storm cell, d(r) centered at the 
bucket. Since the occurrence of storm centers is uniform over the 
universe of storm center locations, the probability that a storm 
center occurs within this circle is proportional to the ratio of the 
area of this circle, A D (r), to the area of an unknown universe, A. 
in which storms are assumed to occur. Thus, the probability that 
a storm center occurs in A g (r) is 
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(26) p(storm in A g ) = A g (r)/A u . 

It follows that the tipping bucket rain rate probability density 
function in terms of the model storm cell rain rate probability 
density function is 

(27) P B (r) = p f (r) A g (r)/A u . 

Due to the physical nature of the precipitation process, the 
function p g (r) has two components. The first, p g (r) is the probability 
density of a non-zero rain rate r in A g (r). The second is p gQ 6(r) 
which is the probability density of no cell center in A g (r). 

Substituting these components for P B (r) in Eq. (27) and 
solving for P r (r) the result is 

(28) p r (r) = [p B (r) + P Bo <s(r)] A u /A g (r). 

The probability distribution function (a Q ,r)) is 
determined from the storm cell propagation path intersection geometry. 
The area in which a storm cell center can occur and intersect the 
propagation path over a length s,^(a Q ,r) is denoted by AU-j (a Q ,r)) ; 
then the probability function p.(ii, (a ,r)) is given by 

I O 

(29) P £ U 1 (a 0 ,r)) = A(£ 1 (a 0 ,r))/A u . 

Substituting Eqs. (28) and (29) into (24) yields 
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CO 


(30) P a (a Q ) = N | [pg(r) + P Bo ( r ) A ^l^ a 0 » r ) )/ A B ( r ) dr ' 

W“> 

The constant N found from Eq. (25) is 


(31) N = 1/ 


P B o A (*i( 0 » 0 ))/ a b(°) + PB^ r ^ A ^l (°> r ))/ A B^ r ^ dr 


The propagation path - storm cell intersection can occur in 
one of four ways. The path can enter the cell side and exit through 
a cell side or through the cell top. On the other hand, the path 
can originate from within the cell and exit through the cell top 
or cell side. Horizontal and vertical views of each of these four 
cases are shown in Fig. 3. 

The boundaries of the area A(^(a Q ,r)) are determined by the 
condition that the intersection of the storm cell and the propagation 
path equals or exceeds a specified length. The restrictions on the 
cell center location (x Q ,y o ) are determined from the equations 
which describe the intersection of the propagation path and the 
storm cell. It is noted first that a necessary condition, which must 
be met, is that the intersection path length specified by £-|(a 0 ,r) be 
shorter than the longest possible intersection path for any cell 
location. The maximum intersection path is shown in Fig. 4 for the 
two exit cases. In Fig. 4a it is seen that the maximum intersection 
length, £ , is determined by the cell diameter, d(r) and the ele- 

ImuA 

vation angle, £. For this case £ max is given by 
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Case a . Side-side cell -path 
intersection 


Case b . Side-top cell -path 
intersection 



Case c . Interior-side cell- 
path intersection 


Case d . Interior-top cell- 
path intersection 


Fig. 3. --Cell-path intersections for a single site. 
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Maximum path intersection length 
determined by the cell diameter. 



h 



Maximum path intersection length 
determined by the cell height. 


Fig. 4. --Maximum path intersection lengths. 
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(32) - ‘max = d(r)/cos 5 ‘ 

The maximum length, , for the top exit case is determined by 
the height, h, and elevation angle, £, as 

( 33 ) = h/sin 5 . ‘ 

Thus, given , the first condition which must be met is 

(34) £-j ^minimum (d(r)/cos h/sin E,). 

Otherwise , A(s,.| (a 0 ,r) ) is zero. 

Assuming this condition has been met, the intersections shown 
earlier in Fig. 3 can be considered. For cases (a) and (b) the 
earth terminal is outside the storm cell so that the propagation 
path enters the storm cell side. Thus, the storm center must lie 
outside the circle given by Eq. (35). 

{35) x o + y o = R2 * 

For these calculations, the rain rate, r, will be assumed to be 
constant and the cell radius will be represented as R to simplify 
the notation. In case (a), the propagation path exits the side 
of the storm. For the intersection points to lie on the circular 
boundary of the cell , the x coordinate of the center, x Q , must be 
located outside the circle given by 
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(36) 


(x Q - h/tan k) 2 + y 2 = R 2 . 

The case (b) intersection shows the propagation path exiting the 

storm cell through the top. The conditions on the storm center 

location for this intersection require that x lie within the 

o 

circles 

(37) (x Q - (h/tan C - ^ cos z)) 2 + y 2 = R 2 
and 

(38) (x Q - h/tan z) 2 + y 2 = R 2 . 

Substituting these conditions on x Q into Eq. (35) the 
restriction that y Q must lie between the lines 

(39) y Q = ± Jr 2 -(£ 2 cos 2 ? )/4 
results. 

The boundaries for the storm center locations for cases (c) 

t 

and (d) are determined in a similar fashion. For these cases, the 
earth terminal must be contained in the storm cell, so the center 
(x 0 ,y 0 ) must lie within the circle 

(40) x 2 + y 2 = R 2 . 

For the propagation path to exit through the cell side, x Q must be 
inside the circle 
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(41) (x Q - £ 1 cos %) 2 + y* = R 2 
and outside the circle 

(42) (x Q - h/tan ?) 2 + y 2 = R 2 . 

The top exit occurs when x Q lies within the circle 

(43) (x Q - h/tan + y 2 = R 2 . 

As before, the limits on y Q are defined by Eq. (39), 

When these conditions are combined the region where a storm 
center may be located and intersect the propagation path is bounded 
by the solid line shown in Fig. 5. The cross-hatched area is the 



Fig. 5. --Area for storm cell -propagation path intersection 
length >_ ^ for a single site. 
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region whose area is A(t.j (r,a o )). The area contained in these 
boundaries is given as 


(44) A(4 1 (o 0 ,r)) = 


hD 


3 D^i cos s' 


tan s 


1 - 


2 2 

J^cos s 


where 


D 2 . -1 
+ sin 


■ 2 :or 2 - 

X 


£^COS s 


D 


(45) 


D = d(r) 


(46) 


a-|(a 0 ,r) = a Q /kr p . 


Recalling Eq. (30) 


(30) 


P (aj = N 

a 0 


r • (a ) 
mi n ' o 


[Pg( r ) + P Bo <s(r)] A(a ] ( a Q ,r))/A B (r)dr 


it is seen that the attenuation probability distribution function has 
been calculated for a single earth terminal which is effected by 
a single model storm cell intersecting the propagation path. In the 
next chapter, the effect of a single model storm cell on two spatially 
separated earth terminals will be considered. 
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CHAPTER III 

CALCULATION OF THE TWO TERMINAL JOINT ATTENUATION 
PROBABILITY DISTRIBUTION 

In the previous chapter the single terminal attenuation 
probability distribution function was calculated using a cylindrical 
storm cell model. This same storm cell model can also be used to 
calculate the two terminal joint attenuation probability distribution 
function due to a single storm cell. If the spatial separation is 
not too large the intersection of the two propagation paths by two 
different cells is not likely and thus a single cell model is all 
that is needed. For this calculation the same storm cell model used 
in Chapter II will be used by assuming that this single storm cell 
affects both terminals simultaneously. 

To calculate the probability that the attenuations observed at 
both terminals simultaneously equal or exceed the attenuation, a Q , 
denoted as P 2 a ( a 0 )> the procedure used earlier to calculate p a (a 0 ) is 
used again. Since the only difference in the calculation is the 
simultaneous propagation path intersection condition, Eq. (21) can be 
used directly to calculate P 2 a ^ a 0 ) if the definition for the path 
intersection probability, P^(^(a Q ,r)), is changed. If the attenu- 
ation observed at both sites equals or exceeds, a Q then the model 
storm cell must intersect both propagation paths in such a manner that 


23 



both of the intersection path lengths equal or exceed £ -|( a 0 » r )* 

With this new definition for (a Q >r)) it 1S seen that P 2 a ( a 0 ) is 

00 

(47> P 2a ( “o> = N { + P Bo S(r):l P 2» U l (a 0' r)> 

A u /A B (r) dr 

where 

(48) P 2j^ £ l^ a o‘ r ^ = ProhaMUty distribution function that 

both propagation paths are intersected 
over a length equal to or exceeding 
t'j (a o ,r ) » simultaneously. 

As before, P 2£^i^ a 0 ,r ^ 1S 9 iven b y 

(49) P 2£^l ^ a o ,r ^ = V*1 ( r ’ a 0 )) /A u 

where (r,a Q )) is the area where the storm center may be located 
such that the storm cell intersects both propagation paths over a 
length equal to or exceeding z^{a Qi r), simultaneously. Since 
A(a 1 (a ,r)) is the area where the storm cell intersects a single 
path over length (a Q ,r), then the overlap area of these single path 
intersection areas is the area, (° 0 »r)). A typical case is 

shown in Fig. 6. Because of the symmetry of the areas, AU-j (a Q ,r)) , 
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y. 



Fig. 6.— Overlap areas for two site probability 
distribution calculation. 

for each terminal it is only necessary to consider the case where 
site one is located at the origin and site two is located in the 
third quadrant as shown in Fig. 6. The calculation of the overlap 
area is straightforward since boundaries of the intersection areas 
are known from the single terminal calculation. Thus, it is necessary 
only to identify the various overlap cases, calculate the boundary 
intersections, and evaluate the area. In calculating the overlap 
area it is assumed that site 2 is separated from site 1 by a distance 
p x in the x direction and in the y direction. These separation 
dimensions are also shown in Fig. 6. 
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Due to the circular shape of the cell model, portions of the 
overlap area are truncated circular segments like the shaded area 
in Fig. 7. To simplify the equations for A 2s,^ £ l^ a o’ r ^ the funct ' ion 
A s ( x i, x 2 >w) will be defined in Eq. (50). 



Fig. 7. —Truncated circular segment. 


(50) A s (x-j, x 2 ,w) = | Jr 2 -(x-w) 2 dx. 

X 1 

In the remaining calculations, the dependence of A s (x^,x 2 ,w) on the 
cell radius, R, will be understood. To further simplify the equations 
the following quantities will be used 
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(51a) 

a 

= h/tan x 

(51b) 

b 

= £-| COS i 

(51c) 

c 

= a + p x 

(51 d) 

d 

= b + p x 

(51 e ) 

e 

CM 

-o 

CM 

■ \ 
II 

(51 f ) 

f 

= e + Py 

(51 g ) 

g 

= p y 

(51h) 

h' 

= h/tan 5 

(51 i) 

a' 

= £.j cos 


Substituting into the boundary equations determined earlier for 
A ( a-| (a Q ,r ) ) in Eqs. (34)- (42) the boundaries for the typical overlap 
shown in Fig. 8 are obtained. Since the relative positions of the 



Fig. 8. --Influence area boundary equations. 
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influence areas will remain the same for all cases to be shown, 
these equations will not be shown on later figures of the overlap 
cases. 

As the placement of the second site is varied in the third 
quadrant eight cases must be considered and these are shown in 
Fig. 9. For each case the intersection points are identified along 
with other pertinent points when space and clarity permit. The 
intersection points were calculated for the eight cases and the 
conditions for that type of intersection determined. The overlap 
areas were then broken into subareas and the subareas calculated. 
The total area is the arithmetic sum of the subareas. The case 
numbers, points used in the subarea calculations, intersection 
restrictions, and subareas are presented in Table 1., To calculate 
P 2 (i, (a 0 ,r)) the intersection points are calculated and the ap- 
propriate case chosen from Table 1. The subarea equations are then 
used to calculate ^Ul (cx 0 »r)). In conclusion, using Eq. (47) and 
the appropriate equations from Table 1 it is possible to calculate 
the joint attenuation distribution for two spatially separated 
terminals . 
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CHAPTER IV 

MODEL STORM CELL PARAMETERS 


The attenuation* probability distributions presented in the 
last two chapters were calculated using a cylindrical storm cell 
model. In this chapter storm cell height, h; storm cell diameter 
as a function of rain rate, d(r); and the tipping bucket rain rate 
probability density function, Pg(r),will be discussed. These 
parameters are used to relate the storm cell model to real world 
storm phenomena. 

Rain Rate Probability Density Function 

The most difficult parameter to determine is the tipping 
bucket rain rate probability density function, p B (r), because of its 
temporal and geographic variability. For the storm cell model to 
be useful in making earth terminal design decisions, the tipping 
bucket rain rate probability characteristics for each earth site must 
be known. Currently, the most commonly available data for calculating 
this function are collected by the National Weather Service as 
mentioned earlier. The rain rates which are measured, however, are 
determined on a clock-hourly basis and represent the average rain 
rate during a sixty minute period beginning at the start of each hour. 
This integration time is too coarse for use here since a rain rate 
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of 100 mn/hr for fifteen minutes and a rain rate of 25 nrn/hr for 
sixty minutes are measured as the same clock-hour rain rate of 
25 mm/hr. The measurement that is needed is the instantaneous rain 
rate. 

Bussey[10] made the first attempt to estimate instantaneous 
rain rates by empirically determining a relationship between clock- 
hour rain rates and instantaneous rain rates. Although his results 
are valid only for Washington, D.C., this approach represents a 
useful initial effort and has been used to estimate instantaneous 
rain rate statistics for other ci ti es [1 1 ] . 

A study of excessive precipitation was made by Dyck and 
Mattice[12] in which an attempt was made to catalog the number of 
occurrences of rain events which produced 2.5 inches or more in 
24 hours or 1.0 inch or more in one hour on a national basis. However 
the data presented is much too coarse to be useful in predicting mi 11 i 
meter wavelength attenuation. Lenhard, Cole and Sissenwine[13] 
presented an algorithm for estimating the rain rates exceeded for 
0.1, 0.5, 1.0 and 2.0 percent of time from the local mean temperature 
and a precipitation index. Although this technique gives four points 
on the rain rate probability curve, further detail is needed. 

A rain rate probability density function, which gives a 
functional form for p B (r) which is geographically determined, has 
been postulated by Rice and Holmberg[14]. The rain rate probability 
distribution function was developed for twelve geographic regions 
in the continental United States for averaging times of one, five 
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and sixty minutes. Although not ideal, the clock-one minute distri- 
bution will be used here since it appears to be the finest scale 
presently available. Since the Rice-Holmberg rain rate probability 
function is used extensively in the attenuation probability function 
calculations it will be briefly explained here. 

This rain rate probability function was derived using precipi 
tation data from forty-nine locations in the United States, as well 
as fourteen other stations throughout the world. From this data 
Rice and Holmberg determined that the rain rate probability distri- 
bution for non-zero rain rates could be modelled by the sum of three 
exponential modes. The total distribution is the sum of these three 
modes as shown in Eq. (52) 

(52) P t (r „J ■ j i exp(- r 0 /B 1it ) 

where P t (r Q ) = Probability distribution that the clock-t 

minute rain rate, r Q , is exceeded 

t = the clock averaging period 

A'l "hh 

i ,t = the mode coefficient for the j region and 

the i mode for the clock-t period 

XL 

B. + = the mode rate for the i xr mode and clock-t 
period. 

Since the Rice-Holmberg function is a probability distribution 
function the density function used in Eqs. (24) and (47) can be 

calculated by differentiation. It is noted that since P^(r Q ) is the 
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probability that a non-zero rain rate r Q is exceeded the resulting 
density function will be for non-zero rain rates only. Thus 

(53 > p B (r o> = i I 1 A i.i /B i,i ex P ( - r o /8 i.n ) - 

It follows from Eq. (3) that the probability density of a zero 
rain rate is 


(54) 




Thus, a form for Pg(r) which is geographically dependent is now 
available. The method in which the local rain rate characteristics 
can be incorporated will now be explained. 

The total rainfall for the observation period used in determin 
ing the A:? . can be calculated by integrating Eq. (53) over all 

1 > t 

non-zero rain rates. It is noted by Rice and Holmberg in Reference 
[14] that the result is independent of the clock averaging time, t. 
Assuming that the rain accumulation for each mode is independent of 
t, then the product A^ ^ ^ is a constant for all values of, t. 

Calculations of this product from data presented in Reference [14] 
shows that this is a reasonable assumption. This relationship of 
the total rain accumulation for each mode can be used to calculate 
the mode coefficients for the clock-one minute density function from 
the clock-hourly data published by the National Weather Service in 
the following way. 
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First, Eq. (52) is fitted, using a minimum mean square error 
algorithm, to the clock-hourly rain rate probability distribution 
calculated from National Weather Service data. Using the values 
for B.j gg and -j obtained from Reference [14], 


B 1 

,60 = 

10 


mm/hr 

B 2 

,60 = 

3 

.0 

mm/hr 

B 3 

,60 = 

1 , 

.0 

mm/hr 

B 1 

,1 = 

35 

.0 

mm/hr 

B 2 

,1 = 

10 

.5 

mm/hr 

B 3 

,1 = 

3 

.5 

mm/hr 


the clock one-minute mode coefficients are calculated from 


(55) 


A i,l = B i ,60 A i ,60 calc /B i,l 


It should also be noted that Rice and Holmberg present a complete 
set of coefficients for twelve geographical regions in the United 
States for clock one-, five-, and sixty-minute rain rates which can be 
used in lieu of the coefficients calculated using the above method. 

The above method, however, incorporates the actual precipitation 
characteristics for a particular site location within a geographic 
region and is preferred for calculating P a (a Q ) and P 2cJ a o^' 

Thus, the Rice -Holmberg rain rate probability function provides a 
mathematically tractable function which permits the calculation of 
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the tipping bucket rain rate probability density function for a 
particular site using locally measured tipping bucket rain rate 
data . 

Cell Diameter - Rain Rate Function 

The remaining storm cell model parameters which must be 
specified are the diameter-rain rate relationship, d(r) and the cell 
height, h. Although d(r) is not directly measurable, two forms have 
been proposed. 

The first, was inferred by Hogg[15] from sun-tracker measure- 
ments and is shown in Fig. 10. Since discrete data points are 
difficult to use in numeric computations an exponential function was 
fitted to Hogg's model using a minimum mean square error criterion 
and is given in Eq. (56). 

(56) d(r) =11.34 exp(-r/33) + 1.25 (km). 

The diameter calculated from Eq. (56) is also shown in Fig. 10. 

A second diameter-rain rate function, recommended by CCIR 
and discussed by Benoit[15], is the logrithmic relationship given 
below. 

(57) d(r) = 41.4-23.5 log 1Q (r) (km). 

This model has two problems. First, for small rain rates the cell 
diameter approaches infinity. Second, the cell diameter becomes 
negative for rain rates exceeding 58 mm/hr. To prevent this from 
occurring in the numeric calculations, the minimum cell diameter 
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RAIN RATE,r( mm / hr ) 


Fig. 10. --Diameter versus rain rate relationship 
from Reference [12]. 
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CELL DIAMETER , d (r) , (Km) 


was clamped at 1.25 km. This value was chosen based on the minimum 
cell diameter calculated 'from Hogg's model. A comparison of the 
exponential model and the modified CCIR model is shown in Fig. 11. 



Fig. 11 .--Comparison of two rain rate diameter 
relationships . 
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Cell Height 


The remaining cell parameter, h, specifies the height of the 
attenuating region within the cell. Although, rainstorms have tops 
tens of thousands of feet in the atmosphere, the liquid precipitation 
is confined to the region below the zero degree isotherm. This 
isotherm represents the altitude at which precipitation is turning 
from solid to liquid form. Several estimates of this height have 
been made and will be discussed here. The first of these was a 
study of thunderstorms in Florida and Ohi o[l 7] . From graphical data , 
presented in this study, an average cell height appears to be approxi- 
mately 6 km. Another source, the CCIR, recommends that a cell height 
of 7 km be used for temperate climates[15], and a paper recently 
presented by Crane[18] on backscatter due to rain used cell heights 
of 15 km in Florida, 11 km in New England and 7 km for the Pacific 
Northwest, California and the desert. It will be seen later from 
calculations of the attenuation probability distribution that the 
results are relatively insensitive to variations in the cell height 
and, hence, this parameter is not considered to be critical. The 
values used in later calculations will be specified for each case. 

Attenuation -Rain Rate Relation 

A final relationship which must be considered is the 
attenuation-rain rate relationship used in Eq. (5) and repeated below. 

(5) a = kr P £ 
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where 

k and p are frequency dependent constants 
jl = path length through precipitation in km 
r = rain rate in mm/hr 
a = observed attenuation in dB. 

Calculations and measurements of k and p can be found in References 
[3], [5], [7], [15], [19], [20] and [21] with a diversity of numerical 
values. The summary of results prior to 1965 presented by Medhurst[19] 
showed an 8 dB variation in the calculated attenuation for a rain rate 
of 100 mm/hr. A comparison of several results by Ippolito[21] 
resulted in his choosing k = .035 and p = 1.155 at 15.3 GHz. Interpo- 
lating data presented by Crane[20] from measurements of 4741 rain drop 
size distributions yields p = 1.10 and k = .035. This value will be 
used for calculations presented in the next chapter. 

Summary of Parameters 

Concluding, then, the parameters which will be used in the 
attenuation probability calculations in the next chapter are 

3 i 

(58) P B (r) = l A 3 j/B i j exp(-r/B.^), 

where A;? , and B. ■, are chosen from Ref. [14] or from 

measured rain rate data if available 

(59) d(r) = 1.250 + 11.340 exp(-r/33) (km) 
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or 


(60) 


(61) 

(62) 


d(r) = 


41.4-23.5 log 1Q (r)(km) 
) 1.250 km 


h = 6000 m to 17000 m 
a = .035 r 1 - 10 a(dB). 


r <^ 52 rnn/hr 
r > 52 mrn/hr 
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CHAPTER V 

COMPARISON OF CALCULATED AND EXPERIMENTAL RESULTS FROM ATS-5 


Attenuation probability distribution functions have been 
calculated in earlier chapters for one and two terminal configurations 
using a cylindrical storm cell model. In this chapter the results of 
measurements made using the Applications Technology Satellite, ATS-5, 
15.3 GHz downlink at two locations will be compared to storm cell 
model calculations incorporating the storm cell parameters presented 
in the previous chapter. For both cases considered instantaneous 
rain rate data was available and was used to calculate the Rice- 
Holmberg distribution coefficients. Also presented in this chapter 
are discussions of the sensitivity of the results to a ± 10 percent 
variation in the numerical value of the parameters and a three mode 
exponential diameter-rain rate function. 

Two ATS-5 receiving stations were considered. The first 
station is located at the Ohio State University in Columbus, Ohio, 
and consists of two terminals separated by 

p = -2.74 km 

M x 

py = -2.99 km 

which yields an intersite spacing of approximately 4.0 km. The 
terminals are discussed in References [22] and [23] and the data 
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acquired in 1970 which is used here is discussed in Reference [24]. 

The instantaneous rain rate data was extracted from recording tipping 
bucket rain gauge charts and the Rice-Holmberg coefficients calculated 
The second station considered was located in Rosman, North Carolina, 
and was operated by NASA Goddard Space Flight Center. The measured 
data were analyzed by the Westinghouse Corporation[25] for the 1970 
data period and consists of attenuation data and instantaneous rain 

j 

rate data taken during attenuation measurements. A Rice-Holmberg 
rain rate distribution was determined from this instantaneous rain 
rate data for use in the attenuation calculations. 

The measured rain rate distribution and the Rice-Holmberg 
model for the Ohio State University data are shown in Fig. 12. The 



Fig. 12. --Rain rate distribution at Ohio State University 
during ATS-5 Data Periods in 1970. 


* 
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measured and calculated attenuation probability distribution 
functions for the 1970, Ohio State University single terminal and 
diversity terminal, ATS-5 measurements are shown in Figs. 13 and 14. 



Fig. 13. --Ohio State University ATS-5 single site measured 
and calculated attenuation distribution. 
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ATTENUATION UB) 


Fig. 14. --Ohio State University 1970 diversity site measured 
and calculated joint distribution. 
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Both the exponential diameter function and the modified CCIR function 
were used in the calculations and each distribution is shown. 

Two additional quantities, the diversity improvement, n(a Q ) 
and the diversity gain, G(a Q ), which are defined in Eqs. (63) and 
(64) are shown in Figs. 15 and 16 respectively. 

(63) n<« 0 > " P a ( “o )/P 2a ( “o ) 

(64) G(o 0 ) = a 0 " a i 

where P (a ) = P, (a,). 

a o za 1 

From these figures it is seen that calculations using the storm cell 
model agree fairly well with experimental measurements. This general 
agreement is also seen in calculations using the NASA Rosman instan- 
taneous rain rate data shown in Fig. 17. The calculated and 
measured attenuation probability distributions are shown in Fig. 18. 

The agreement between the calculated and measured attenuation 
distributions is summarized in Fig. 19 for the exponential diameter 
function and the modified CCIR diameter function in Fig. 20. The error 
shown is the difference between the calculated attenuation and the 
measured attenuation for a given value of p a (« 0 )- Thus, this error is 
the horizontal distance between the experimental and calculated curves. 
It is seen in this figure that the error curves have a similar nature 
and range from -3 dB to +3 dB. The magnitude of the error is encourag- 
ing and shows that the storm cell model technique for predicting 
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DIVERSITY IMPROVEMENT 7 ] la) 



SINGLE TERMINAL ATTENUATION (dB) 

Fig. 15. --Diversity improvement for Ohio State 
University ATS-5 1970 diversity data. 
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DIVERSITY GAIN ( d B ) 



SINGLE TERMINAL ATTENUATION ( dB ) 


Fig. 16. --Diversity gain for Ohio State University 
ATS-5 1970 diversity data. 
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MINUTES ATTENUATION WAS EQUALLED OR EXCEEDED 



Fig. 18. --NASA ATS-5 1970 measured and calculated 
attenuation distribution. 
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Fig. 19. --Comparison of calculated and measured 15.3 GHz 
ATS-5 1970 data using the exponential diameter 
function. 
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Fig. 20. --Comparison of calculated and measured 15.6 GHz 
ATS-5 1970 data using the modified CCIR 
diameter function. 


mil limeter wavelength attenuation statistics is potentially useful. The 
similarity of the error curves suggests that the present model has a 
systematic error that might be reduced by modifying the model parameters. 

To examine the sensitivity of the model to variations of the 
model storm cell parameters a ten percent perturbation was imposed 
and the single site yearly attenuation probability distribution was 
calculated for the site located in Columbus, Ohio, using the Rice- 
Holmberg coefficients for the rain rate distribution function. It 
was found that there was little change in the calculated distribution 
due to elevation angle or the cell height changes. However, pertur- 
bations of the remaining parameters, d(r), k, and p had a significant 
effect on the probability distribution as shown in Figs. 21-25. 

The parameter which showed the greatest sensitivity was p in 
Eq. (5), the Gunn-East attenuation relation. The ±10% perturbation 
caused the calculated value of P(a Q ) to change more than 3 orders of 
magnitude for large attenuations. For a constant P(<* 0 ), this 
corresponded to a ±5 dB change in a Q . The high sensitivity of the 
attenuation distribution to variations in p is probably due to 
several factors. First, p is an exponential constant greater than 
one, thus, the perturbations are exponentially magnified. Second, 
p influences both the lower limit of the r integration, r^ n (a), 
and the minimum intersected path length, £^(a Q ,r). 

The perturbations of the other constant in the Gunn-East re- 
lation, k, and the constants in the exponential cell diameter-rain 
rate function, d(r), were less severe. As seen in Figs. 22-25 the 
±10% perturbation of the above named parameters resulted in an order 
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REFERENCE 

O + 10% 

A -10% 

REFERENCE PARAMETER 
VALUES 

h = 6 km 
£ = 40° 

d (r ) = 11.34 exp (~ r / 3 3 ) + 
a (r ) = 0.035 j? r '•'<> 







Fig. 25. --Effect of ±10 percent variation of 
cell diameter rate. 
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of magnitude variation in P(a Q ) and at constant P(a Q ), about a 
±2 dB variation in a Q . Although these variations are significant, 
it is clear that the most critical constant is p in Eq. (5). 

Ultimately the factor which will limit the accuracy of this 
model for use in terminal design is the yearly rainfall variation. 

To examine the magnitude of this variation, four years of rainfall 
data from the Port Columbus station of the National Weather 
Service were used to calculate the annual Rice-Holmberg clock-minute 
rain rate probability density functions. These coefficients were 
then used to calculate the attenuation probability distribution 
functions which are shown in Fig. 26. For reference the total 
yearly accumulations are also shown for each year. The largest 
variation of a Q at a constant P a (a 0 ) is approximately ±1.5 dB. 

While this test does not necessarily represent a worst case it does 
provide a feeling for the order of magnitude of the variations that 
might be expected. 

To improve the accuracy of the calculated attenuation 
probability function at the high attenuations, modifications to the 
model parameters were considered. The parameter which has the least 
certainty associated with its value is the diameter rain rate function 
and this function was examined. Following the idea presented by 
Rice and Holmberg, that three exponential modes were needed to ap- 
proximate the rain rate distribution function, a three mode exponential 
function was postulated for the diameter rain rate function. Using 
the exponential diameter function as a starting point, the calculated 
attenuation probability distribution function for the Ohio State 
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ATTENUATION ( dB ) 


Fig. 26. --Effect of yearly weather variations for 
1965, 1967, 1968 and 1970 at Columbus, 

Ohio. 

University single terminal was fitted by trial -an d-error to the 
measured data by modifying the diameter function parameters. The 
measurements above 8 dB attenuation were ignored since they were 

dominated by the system threshold. The resulting three mode diameter 
function is given in Eq. (65). 

( 65 ) d(r) = 1134. exp(-r/3.5) + 1.13 exp(-r/33) 

+ 1.25 exp(-r/264) (km). 

For comparison the diameter function in Eq. (65) is shown in Fig. 27 
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Fig. 27. --Comparison of three rain rate diameter relationships. 

with the diameter functions mentioned earlier. It is seen that the 
three mode diameter function is similar to the CCIR relation with two 
exceptions. First, Eq. (65) is approaching a finite but large value 
at r=0. Second, Eq. (65) decays exponentially for high rain rates. 

The resulting attenuation probability distributions are shown 
in Fig. 28. It is seen that the general shape of the single terminal 
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Fig. 28. --Ohio State University ATS-5 single site measured and calcu- 
lated joint distribution for modified diameter function. 


probability distribution has been approximated if the threshold 
flattening of the experimental data for large attenuations is ignored. 

Using this new diameter function the diversity attenuation 
distribution was calculated and is shown with the measured distri- 
bution in Fig. 29. It is seen that good agreement has also been 
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Fig. 29. — Ohio State University ATS-5 diversity site measured and 
calculated joint distribution for modified diameter 
function. 


obtained here. As a further test of the new diameter function, the 
attenuation distribution was calculated for the NASA terminal and 
the results shown in Fig. 30. The agreement between the calculated 
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ATTENUATION (dB) 


NASA ATS-5 measured and calculated attenuation 
distribution for modified diameter function. 






and experimental distributions has been improved here also. The 
errors are summarized by the error curves in Fig. 31. It is seen that 


o 

h 

(/> 



Fig. 31 .--Attenuation error for three mode 
diameter function. 


using the three mode diameter function the range of the error in 
predicting the attenuation is between +1 and -2 dB. Therefore, milli- 
meter wavelength attenuation distributions for single terminal and 
two terminal diversity configurations can be calculated with reasonable 
accuracy using this simple storm cell model. 
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CHAPTER VI 
SUMMARY 


In conclusion, a cylindrical, single storm cell model was 
postulated. The attenuation probability distribution for an earth - 
space propagation path was calculated for both a single earth terminal 
and a spatial diversity terminal configuration. The calculation 
procedure used a tipping bucket rain rate distribution function, a 
storm cell diameter- rain rate function, a fixed cell height and the 
Gunn-East attenuation-rain rate relation as parameters. The tipping 
bucket rain rate distribution function proposed by Rice and Holmberg 
was discussed and its usefulness is making calculations for a par- 
ticular geographic location from tipping bucket rain rate measure- 
ments was discussed. A cell diameter-rain rate function measured by 
Hogg and a diameter function proposed by CCIR were presented. Also 
discussed were typical values for the model storm cell height. 

The attenuation distribution functions for two ATS-5 Millimeter 
Wavelength Experiment receiving stations were calculated using the 
chosen cell parameters and measured tipping bucket rain rate data. 

The calculated distributions agreed within 3 dB of the measured 
distributions. The sensitivity of the calculated results to parameter 
variations was discussed and it was determined that the rain rate 
exponent in the Gunn-East relation was the most critical parameter, 
causing a 5 dB variation in the prediction for a 10% variation. 
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It was also found that yearly precipitation variations could cause at 
least a 1.5 dB variation in the attenuation distribution. 

To improve the accuracy of the calculated attenuation distri- 
bution, a three mode exponential diameter function was postulated 
and the Ohio State University single site measured attenuation 
probability distribution was used to determine the diameter function 
coefficients. The three mode diameter function was then used to 
calculate the Ohio State University diversity attenuation distribution 
and the NASA Rosman single site attenuation distribution. The cal- 
culated distributions agreed within -2 to +1 dB of the measured 
distributions which represented a significant improvement over the 
other diameter functions used. 

Thus, it has been shown that a simple cylindrical storm cell 
model is useful in predicting millimeter wavelength precipitation 
attenuation. The model proposed is mathematically tractable and 
sufficiently general so that modification can be easily made. 
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